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ABSTRACT 



We present the redshift lower limit oi z > 0.6035 for the very-high-energy 

W : (VHE; E >100 GeV) emitting blazar PKS 1424+240 (PG 1424+240). This limit 

is inferred from Lyman /3 and 7 absorption observed in the far-ultraviolet spectra 

from the Hubble Space Telescope/Cosmic Origins Spectrograph. No VHE-detected 

blazar has shown solid spectroscopic evidence of being more distant. At this dis- 



O. 

^ . tance, VHE observations by VERITAS are shown to sample historically large 

j^ . gamma-ray opacity values at 500 GeV, extending beyond r=4 for low-level mod- 
els of the extragalactic background light (EBL) and beyond r =5 for high-levels. 

CN , The majority of the z = 0.6035 absorption-corrected VHE spectrum appears to 

Q\^ . exhibit a lower flux than an extrapolation of the contemporaneous LAT power- 

^ . law fit beyond 100 GeV. However, the highest energy VERITAS point is the only 

Tj- , point showing agreement with this extrapolation, possibly implying the overes- 

■^ . timation of the gamma-ray opacity or the onset of an unexpected VHE spectral 
O 
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feature. A curved log parabola is favored when fitting the full range of gamma- 
ray data (0.5 to 500 GeV). While fitting the absorption-corrected VHE data 
alone results in a harder differential power law than that from the full range, 
the indices derived using three EBL models are consistent with the physically 
motivated limit set by Fermi acceleration processes. 

Subject headings: BL Lacertae objects: individual: PKS 1424+240 (PG 1424+240), 
galaxies: active, intergalactic medium, ultraviolet: general, cosmology: diffuse 
radiation 



1. Introduction 



Blazars are active galaxies with a jet closely aligned with the Earth line of sight 
( lUrry fc Padovanil Il995l ). These objects are the most commonly detected type of extra- 
galactic source at very high energies (VHE; E >100 GeV). The VHE gamma rays which 
propagate through the intergalactic r nedium are absorbed by extragalactic ba ckground light 



(EBL) photons via pair production (JNikishovl Il962l : iGould Sz Shreder 



difficult to directly measure due to strong foreground sources fJHauser fc Dwekl 120011 ) 



19671). The EBL is 



Various methods are used to estimate the density of the EBL. A selection of EBL models 
include a semi-analytical model ( iGilmore et al.ll2012[ ). an observationally based model which 
utilizes observations of i^-band rest frame ga laxy luminosity functio ns in combination with 
galaxy spectral energy distribution fractions (IDominguez et al.ll201ll) and a model based on 
starlight emission and dust re-emission templates ( jFinke et al.l 120101 ) . Each of these est i- 
mate photon densities are above the lower limits set by galaxy counts (IWerner et al. Il2004r). 



and b elow the limits set by detection of extragalactic VHE photons (e.g., lAharonian et al. 



The interaction between VHE and EBL photons produces a "gamma-ray horizon," lim- 
iting the distance to which VHE sources can be detected. 3C 279 is the most distant VHE 
emitting blazar with a spectroscopically measured redshift. By spectroscopically measured, 
we mean the detection of emission and/or absorpt ion lines of the host galaxy. 3C 279 is a 
fiat-spectrum radio quasar (FSRQ) at z = 0.536 (IMarziani et al.l 119961 ) and was observed 
above 50 GeV by MAGIC during a flare jAlbert et al.lboosi ). 

FSRQs are blazars which show emission lines with an equivalent width of > 5A , pro- 
viding strong host-galaxy emission and absorption for spectroscopic distance measurements. 
Conversely, BL Lacertae-type blazars display weak or no e mission/absorption li nes, and lack 



evidence for a Ca H+K break in their optical spectrum fJMarcha et al. 



19961 : iHealev et al. 
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20071 ). With this definition, it is not surprising that there are VHE-detected BL Lac objects 
that lack definite spectroscopic redshift measurements. Some notable VHE blazars with- 
out spectroscopic redshift are KUV 00311-1938, PG 1553+113, 3C66A, S50716+714 and 
PKS 1424+240 (PG 1424+240). 



KUV 00311-1938, detected at TeV energies by HESS (IBecherini et al. 1120121 ). has a 
tentative redshif t of 0.61 from the observa tion of five features in a poor signal-to-noise ratio 
(SN) spectrum (jPiranomonte et al.l 120071 ). Follow-up spectroscopy with improved SN did 
not confirm any z = 0.61 features, although an intervening Mg II absorption system at 
z = 0.506, not v isible in the prev ious spectrum, was observed and sets a lower limit on the 
blazar distance ( Pita et al.ll2012l ). 



Since BL Lac objects display bright and featureless spectra, strict lower limits on the 
redshifts can be inferred through observation of far-UV (FUV) absorption by the low-z 
intergalactic medium (IGM). This method has been applied to PG 1553+113, 3C66A, an d 
S5 0716+714 s howing lower limits of z =0.3 95, 0.3347 and 0.2314 JPanforth et al.l J2010h : 



Furniss et al.l (120131 ): iDanforth et al.l (l2013al ). respectively). We apply this technique to the 
VHE-detected blazar PKS 1424+240, constraining the blazar to 2; > 0.6035, making it the 
most distant VHE-detected blazar thus far. 



2. Observations and Far-UV Spectral Analysis 

PKS 1424+240 was observed under a Hubble Space Telescope (HST) program (12612, 
PI: Stocke) which uses fiaring blazars as probes of intervening, weak IGM absorption. 
PKS 1424+240 is one of ~ 200 objects monitored b y a network of automate d telescopes: 
the Katzman Automatic Imaging Telescope (KAIT, iFilippenko et al. I (1200 ll )). the "NF/ 



Observatory" (IGrauer et al.ll2008f) a nd th e Small and Moderate- Aperture Remote Telescope 
System fSMARTS. iBonning et all J2012[ )l Optical photometry in April of 2012 indicated 
that PKS 1424+240 was sufficiently bright to trigger a five-orbit HST/COS observation. 

The blazar was observed on 19 April 2012, with the medium resolution (Af ^ 18 km s~^), 
far-UV gratings G130M (1135 < A < 1450 A, 6.4 ksec) and G160M (1400 < A < 1795 A, 
7.9 ksec). The fiux-calibrated, one-dimensional spectra were obtained from the Mikulski 
Archive for Space Tel escopes and combined with the standard IDL procedure s descri bed in 
Danforth et al.l (120 lOl ). Additional analysis details are given in lDanforth et al.l (J2013bl ). The 
combined data show a continuum fiux level of ~ 1.4 x 10~^^ erg cm~^ s~^ A~^ and a median 
SN per 7-pixel resolution element of ~ 19 over the entire spectrum. The spectrum was nor- 
malized with an iterative, spline-based procedure and an automated line-finding algorithm 
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was used to locate > 4a absorption features (JDanforth et al.ll2013bf ) 



We set a firm, observational lower limit to the source redshift by examining intervening 
absorbers. We see absorption consistent with higher-order H I Lyman absorption at z ^ 0.6 
(Figure 1). The presence of absorption profiles consistent with Ly/3 and Ly7 at three distinct 
redshifts provides unambiguous line identifications of absorbing gas along the line of sight 
(Table 1). The measurements summarized in Table 1 are made using Voigt-profile fits to 
the lines. The significance levels are calculated from the o bserved equiv a lent w idths and 
flux errors and are calculated via the methods described in iKeeney et al.l (120121 ). The line 
identifications are bolstered in two cases by low-significance Lj6 detections. 

In other cases where H I absorption is used to constrain blazar distance, the lack of 
absorbers past a certain redshift is used to place a statistical upper limit. For PKS 1424+240, 
the Lya forest extends to the red end of the FUV detector (~ 1800 A, z = 0.47). Absorbers 
at higher redshift are detected through paired Lj(3 and Ly7 lines visible at z < 0.75. However 
this technique is less sensitive to H I absorbers than detection via Lya since both Ly/3 and 
Ly7 lines must be detected to unambiguously identify an absorber and (/ A)Ly7 is only ~ 5% 
that of Lya. 

No H I absorbers are seen between the reddest system [z = 0.6035) and the red edge 
of the detector in Ly/3 [z = 0.75) for a line-free path length Az = 0.15. The spectral data 
quahty is such that we should detect lines at log A^h i ^ 14.0 (VFLy^ ^ 60 mA, Wi^yj > 20mA) 
at a 4(j level. The freque ncy of lines of this strength or higher at low redshift is dH /dz ^ 24 
(JDanforth fc Shullll2008l ). so we would expect M = 3.5^^;8 absorbers to be present if the 
source were at 2; > 0.75. Since no lines are seen, we can rule out z > 0.75 at approximately 
2cr confidence. More detailed simulations may refine the redshift upper limit, but near-UV 
spectra and a direct search for Lya lines at 2; ~ 1 will be less ambiguous. 

The redshift low er limit is signif icantly higher than the co nservative estimates of z > 0.06 
( Scarpa et al.lll995l ) and z = 0.23 ( Meisner fc Romani I2OIO), der i ved fr om host-galaxy as- 
sumptions. Notably, the redshift estimates in iMeisner fc Romanil (J2010l ) are in close agree- 
ment with other blazar distance limits derived from the observation of absorption systems. 
The redshift limit for PKS 1424+240 is also below the upper limits oi z = 0.66 and 1.19 
derived from correcting VHE observations f or EB L ab sorption to ma t ch th e Fermi Large 
Area Telescope (LAT) data in lAcciari et al.l (I2OIOI ) and lYang fc Wand (120 lOl ). respectively. 
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3. Absorption of Very High Energy Photons 

The energy- and redshift-dependent absorption of VHE gamma rays by the EBL can 
be estimated using the model-specific optical depth, t{E, z), where the absorption-corrected 
(deabsorbed) flux, F^or, is estimated from the observed flux, Fof,s, using the relation Fcor = 
Fobs X e'^'^^'^\ We investigate the deabsorbed VHE spectrum of PKS 1424-1-240 at the redshift 
lower limit oi z = 0.6035, using three models to explore the effect of relatively low, medium, 
and h igh levels of EBL photon density. The lowest density model used is from lGilmore et al. 
( I2OI2I ). which estimates a. z ^ EBL spectral energy dis tribution nearing the req uired lower 



limits on the EBL set by galaxy counts. We use the iDominguez et a 



estimate VHE absorption by an intermediate level EBL density and the 
model to represent a relatively high level of EBL density. 



301 



(20111) model to 



Finke et all (120101 ) 



3.1. Constraining the Opacity of the EBL 

The blazar VHE spectral index can be used to estimate the EBL spectral properties 
under the assumption that the intrinsic spectrum, cha racterized by the power- law dN/dE ex 
E~^ , cannot be harder than F = 1.5, as described in lAharonian et al.l ( 120061 ). The limit is 
physically motivated by the shock acceleration paradigm, where the hardest index obtained 
for the accelerated leptons is 1.5. 

Notably, the hardest blazar spectral index measure d by the L AT has an index of 1 . 1 but 
is in statistical agreement with the theoretical limit ( lAckermann et al.l 1201 It iNolan et al. 
2OI2I ). Since the LAT is most sensitive to photons at energies where EBL absorption is 
negligible, the indices derived from LAT observations reflect the intrinsically emitted gamma- 
ray spectrum. A more conservative limit equal to the LAT measured index can be placed 
under the assumption that blazars do not harden with inc reasing energy . For PKS 1424-1-240, 
the contemporaneous LAT-measured index is 1.80±0.07 (jAcciari et al.ll2010l ) (Figure 3). 



A power-law flt can be applied to the absorption-corrected points for the redshift lower 
limit oi z = 0.6035, as summarized in Table 2 for each of the EBL models. The flt- 
ted indices for the d eabso rbed spectra usin g the relatively low and medium EBL models 
( JGilmore et al.l (120121 ) and IDominguez et al.l ( 1201 ll ). respectively) are well within the F=1.5 
and F = 1.80±0.07 limits. However, the F = 0.6±0 . 8 inde x resulting from deabsorption with 



the relatively more opaque model from lFinke et al.l ( I2OIOI ) is below, but still consistent with, 
either of the expectations. Improved gamma-ray observations may reveal that this model 
is too dense. Other explanations, such as time-depende nt stochastic accelerated inverse- 
Compton scenarios ( iLefa et al. 1 1201 lal : iLefa et al.l l2011bf ) and internal gamma-gamma ab- 
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sorption (jAharonian et al.ll2008l ). can also account for unusually hard VHE emission spectra. 



3.2. The Gamma-ray Horizon 

Intergalactic gamma-ray opacity due to the EBL has direct consequences for the es- 
timation of the intrinsic gamma-ray spectra of extragalactic VHE targets, with the source 
emitted flux being suppressed by e~^(^'^\ This energy- and redshift-dependent flux suppres- 
sion requires sources to be exponentially brighter at larger distances in order to be detectable 
at VHE. 

The opacities probed through the VHE observation of blazars with redshift informa- 
tion provides insight into the possibility of a pair-production anomaly, as investigated in 
Horns &: Meyerl (120121 ). In the seven VHE blazar spectra which probe opacities in the range 
1< r < 2, an upturn of the absorption-corrected spectra is apparent with a significance 
of 4.2(7 at the r > 2 transition energy. Due to the different energies of the r > 2 transi- 
tion for the blazars studied, source-intrinsic features are an unlikely explanation. This study 
was limited by the number of VHE blazars probing opacities r > 2 with known redshift. 
PKS 1424-1-240 can now be included in the study, expanding the limited opacity parameter 
space available. 

Before the determination of the redshift lower limit of PKS 1424-1-240, the highest sam- 
pled gamma-ray optical depth probed was associated with the detection of 3C 279 during an 
ele vated state, probing a t(E = 475 GeV, z = 0.536) = 3.2, 3.9 or 4.3 when estimated with 



the lGilmore et al.l ( 12012| ): iDominguez et al.l (120111 ): JFinke et al.l ( 120101 ) models, respectively. 
Now, wit h PKS 1424-1-240, th e hig hest opacity s ample d is between 4.1 and 5.3 (if estimated 
with the JGilmore et al.l ( 12012! ) and lFinke et al.l (120101 ) models, respectively.) 



We show the highest opacity sampled by the VHE detection of 3C279 compared to 
the opacity probed wit h the detection of PKS 1424-1-240 in Figure 2, as derived for the 
Dommguez et al.l (1201 ll ) model. The spectral points are plotted along with the t{E, z) = 1—5 
horizons. For reference, we also show published spectral points of every VHE-detected blazar 
with a spectroscopically measured redshift above 0.2 (four sources, as compared to the more 
than 30 VHE blazars with z < 0.2). Since the redshift of PKS 1424-1-240 represents a lower 
limit, the maximum opacity sampled by the VHE detection is illustrated with a right-pointing 
arrow. 
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4. Absorption-Corrected Gamma-ray Emission 



The observed gamma-ray peak is reproduced from lAcciari et al.l (120101 ) in Figure 3, 
illustrating the contemporaneous LAT and VERITAS data with black squares and circles, 
respectively. Additionally, we correct the observed VHE spectrun a for EBL absorption a t 
the minimum redshift oi z = 0.6035 using the intermediate model (JDominguez et al.ll201ll ). 
illustrated by grey circular points. Absorption of the LAT data is negligible. 

The LAT and VERITAS gamma-ray observations extend from 0.5 to 500 GeV. This 
full range, when corrected for EBL absorption sX z = 0.6035, is not well fit with a power 
law (dashed blue line), while a curved log parabola provides an improved fit (dotted blue 
line). We note than neither the power-law nor the log-parabolic fit represents the data well. 
The fit results for each EBL model are summarized in Table 2. Both the power-law and log- 
parabolic fits find indices of F ~ 2, with a significantly curved log-parabolic fit. Notably, the 
highest energy VERITAS point at 500 GeV appears to deviate from the fits, matching the 
LAT extrapolated power-law (red dashed line). This deviation is not expected by standard 
blazar emission models, but is only about two standard deviations at the redshift lower limit. 

A break between the LAT and VERITAS absorption-corrected data is apparent. This 
discrepancy is not likely an issue of instrumental cross-calibration, as agreement between 
VERITAS a nd LAT observations ha ve been found for other contemporaneous blazar obser- 
vations (e.g., lAliu et al.l ( 120111 . |2012[ )). A portion of this feature may be due to a small level 



of undetected variability. Althou gh short intervals of variability are difficult to rule out, no 
long-term variability is detected (jAcciari et al.ll2010l ). making it unlikely that the spectral 
feature between the LAT and VERITAS instruments is due to variability alone. 

Since z = 0.6035 is a lower limit, it is conceivable that the discontinuity between the 
LAT and VERITAS data may, in fact, be an unphysical effect arising from the incomplete 
correction for absorption by the EBL. The first differential flux point of the VHE spectrum 
at 150 GeV cannot be made to match the LAT extrapolated spectrum without deabsorbing 
the flux for a redshift of 1.2 (blue stars in Figure 3). This deabsorbed spectrum is shown for 
the iDominguez et al.l (120111 ) model, b ut is representa t ive of the r equired distances o f z = 1.5 
and z = 1.0 when corrected with the iGilmore et al.l (120121 ) and iFinke et al.l (120101 ) models, 
respectively. 

The 2; = 1.2 corrected VHE spectrum results in a rising slope, i.e. with an index 
F = —2.5 ±1.0 when fit with a differential power law. A VHE spectrum with a rising 
power law is difficult to produce in even the most extreme emission sc enarios. Although thi s 
redshift value is still in agreement with the redshift upper limit set by lYang &: Wang] (120101 ). 
we interpret the unphysical VHE spectrum as evidence that the blazar does not reside at 



this distance. 



4.1. Possible Signature of Intrinsic Gamma-ray Absorption 



Assuming that the blazar resides near z = 0.6035, the apparent discontinuity be- 
tween the LAT and VERITAS energy ranges may be due to gamma-ray absorption in 
the vicinity of the blazar. It has been shown that absorption of gamma-rays by a broad 
l ine region (BLR) can p roduce broken power-law spectra in bright LAT-detected blazars 
(jPoutanen fc Sternll2010l ). However, this type of absorption is not immediately expected for 
an intermediate/high-synchrotr on-peaked source such as PKS 1424-1-240, expect ed to exhibit 
a clean radiation environment (JBottcher &: Dermerll2002l : iGhisellini et al.lll998l ). 



The absorption-corrected VHE point at 500 GeV matches the LAT power-law extrapo- 
lation, with a distinct mismatch to the 100 through 400 GeV points. A source of gamma-ray 
opacity that is only sensitive to photons between 100 GeV and ~400 GeV is difficult to 
explain with an ion continuum such as that present in a BLR. It has been shown that 
the optical depth of a BLR containing UV-continuum an d ionization lines produce s a con- 
stant optical depth from tens of GeV to beyond 30 TeV ( Tavecchio fc Mazinll2009l ). Since 
PKS 1424-1-240 is not expected to harbor a BLR, it is perhaps more likely that the EBL 
model is slightly overcorrecting for the photon absorption around 500 GeV. The observed 
hard spectrum might also be expla ined by secondary garnma rays produced in c osmic-ray 
interactions along the line of sight (JEssey fc Kusenkdl2010l : lAharonian et al.l 120131 ). 



Conclusions 



We present the strict redshift lower limit oi z > 0.6035 for PKS 1424-1-240, set by the 
detection of Ly/3 and 7 lines from intervening hydrogen clouds. This lower limit makes 
PKS 1424-1-240 the most distant VHE-detected source. At this distance, VHE observations 
of the source out to energies of 500 GeV probe gamma-ray opacities of up to r ~5. 

An investigation of possible constraint on the opacity of the EBL shows that the 
absorption-corrected power-law fits do not lie sig nificantly ou t side o f the standard spec- 
tral limitations. However, deabsorption with the iFinke et al.l (J2010l ) model produces the 
hardest intrinsic VHE spectrum. If the blazar resides at a redshift beyond the lower limit, 
the deabsorbed indices may become constraining to even the lowest level EBL model. 

This redshift information allows the investigation of the EBL absorption-corrected gamma- 
ray emission. Correcting the VHE spectrum for the minimum redshift oi z = 0.6035 shows 
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an unexpected spectral shape. The elevated flux around 500 GeV, while not high signifi- 
cance, may be due to an ov ercorrection of absorp tion by the EBL model. This feature is 



still present when using the iGilmore et al.l ( 12012| ) model (open circles in Figure 3), which 



predicts a low EBL density and therefore might be evidence for a pair-production anomaly, 
such as might occur if VHE photons mixed with axion-like particles. 

The spectral feature occurring between the LAT and VERITAS energy bands cannot be 
reasonably removed by correcting for additional absorption due to a larger distance. Instead, 
deabsorption for a higher redshift to match the LAT and VERITAS fluxes results in a non- 
physical spectrum and an extreme distance of 2; = 1.2. It is possible, although unlikely, that 
the EBL absorption-corrected gamma-ray spectrum of PKS 1424+240 is exhibiting gamma- 
ray absorption within an intrinsic photon fleld such as a BLR. As the most distant VHE BL 
Lac object, probing historically high values of gamma-ray opacity, this source requires addi- 
tional gamma ray observations and a tighter constraint on the redshift to further investigate 
the intrinsic emission. 
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Table 1. Intervening Lyman Absorption Lines 



Line Aoba (A) 


Wr (mA)'^ 


b (kms-i)t' 


log N {cm-2)c 


S.L. (<T)d 


z = 0.5838 system 










Ly/3 1624.6 


89 ±29 


47 ±7 


14.18 ±0.12 


6 


Ly7 1540.2 


73 ±5 


55 ±8 


14.52 ±0.04 


5 


z = 0.5960 system 










Ly/3 1637.1 


180 ±1 


44 ±3 


14.52 ±0.02 


11 


Ly7 1552.1 


60 ±33 


37 ±5 


14.44 ±0.19 


5 


LyS 1515.7 


19 ±6 


31 ±5 


14.3 ±0.2 


~ 2 


z = 0.6035 system 










Ly/3 1644.8 


70 ±11 


14 ±2 


14.15 ±0.06 


7 


Ly7 1559.5 


36 ±10 


15 ±3 


14.25 ±0.08 


4 


Ly<5 1522.9 


~13 


~10 


14.1 ±0.3 


'v. 2 



^Rest-frame equivalent width. 

'^Doppler parameter. 

'^Corresponding column density. 

'^ Line significance (according to lKeenev at al.l 1120121) ). 



Table 2. Result s for the power-law fit to the absorption-corrected VHE points from 
Acciari et al.l ( l2010l ) for z = 0.6035. Additionally, we show the fits to the full range of data 



(0.5-500 GeV) for a power law jdN/ dE oc {E /En) ^) a nd lo g parabola 

{dN/dE oc (E/Eo)-^-^L°g(^)). The fits for the IPommguez et al.l f l201l[ ) model are shown in 

Figure 3 in blue dashed and dotted lines, respectively. 



EBL 


VHE Range 


X^ 


Full Range 


X^ 


Full Range 


Full Range 


X 


Model 


Deabsorbed 


(4 DOF) 


Deabsorbed 


(9 DOF) 


Deabsorbed 


Deabsorbed 


(8 DOF) 


Used 


Power-law 




Power- law 




Log-parabola 


Log-parabola 






Index 




Index 




Index 


Curvature 






r 




r 




r 


/3 





Gilmore et al. f2012) 


1.5±0.8 


0.85 


2.07±0.03 


2.3 


2.04±0.04 


0.10±0.03 


1.4 


Domineuez et al. ("2011) 


1.0±0.7 


0.88 


2.01±0.03 


2.0 


1.99±0.04 


0.08±0.03 


1.5 


Finke et al. f2010) 


0.6±0.8 


0.83 


1.97±0.04 


1.9 


1.96±0.04 


0.07±0.03 


1.5 
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Fig. 1. — COS spectra show intervening absorption systems at z = 0.5838, z = 0.5960, and 
z = 0.6035 (arrows) toward PKS 1424+240 in Ly/3 (1025. 72A), Ly7 (972.54 A), and Lj6 
(949.74 A). COS flux and error (grey) vectors are binned by four pixels (lialf a resolution 
element). Continuum fit is shown with dashed line. Other intervening absorption is iden- 
tified with species and redshift in red. Lya absorption systems are observed to the edge 
of the detector (I8OOA, corresponding to ziya < 0.47). Galactic {v ^ 0) absorption and 
instrumental features are l abeled in green. Lya, and 7 features are detected at > 3a and 
will be further discussed in iDanforth et al.l (J2013bl ). 
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PKS 1424+240 from Acciari et al. 2010 
3C 279 from Albert et al. 2008 
PKS 1 222+1 22 from Aleksic et al. 201 1 
1 ES 041 4+287 from Aliuetal. 201 2 
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x=3 from Dominguez et al. 2011 
T=4 from Dominguez et al. 2011 
x=5 from Dominguez et al. 2011 
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Fig. 2. — The highest energy points of the VHE-detected blazars with pubhshed VHE data 
and spectroscopic redshifts beyond 0.2. The close proximity of VHE blazars is a result of 
the gamma-ray opacity of the Universe. The r = 1 — 5 gamma-ray horizon contours are 
shown as bands, including model errors, representing the energy and redsh ift dependent e" 



suppr ession of the VHE flux for extragalactic sources as calculated from iDommguez et al. 



( I2OIII ). The VHE detection of PKS 1424+240 is shown with a rightward arrow, indicating 
the redshift is a lower limit. 
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Fig. 3. — The gamma-ray peak of the spectral energy distribution of PKS 1424+240, with 
LAT (squares and power - law fit contour) and VERITAS observations (black circles) taken 
from I Acciari et al.l ( l2010l ). An upper limit at 750 GeV is shown with a downward pointing 
arrow. The LAT data have been selected to be contemporaneous with the VERITAS observa- 
tions. T he absorption-corr e cted V HE spectrum is shown with the grey circles, using opacities 
from the iDominguez et al.l (120111) EBL model. For reference, the absorption-corrected points 
using the lGilmore et al.l fl2012l) model are shown in open circles, with errors (not drawn) sim- 



ilar to those shown for the iDominguez et al.l (120111 ) deabsorbed points. The LAT power-law 
fit has been extrapolated up to VHE (dashed red line). Power-law and log-parabolic fits to 
the full range (0.5-500 GeV) are shown in the blue dashed and dotted lines, respectively, 
with fitting results in Table 2. To bring the first absorption-corrected VERITAS spectral 
point to match the LAT observed spectrum, the blazar needs to be corrected for absorption 
expected for z ~ 1.2, shown by blue stars (the upper limit for this deabsorption is off-scale). 



